Local density of states of electron-crystal phases in graphene in the quantum Hall 

regime 
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We calculate, within a self-consistent Hartree-Fock approximation, the local density of states 
for different electron crystals in graphene subject to a strong magnetic field. We investigate both 
the Wigner crystal and bubble crystals with Me electrons per lattice site. The total density of 
states consists of several pronounced peaks, the number of which in the negative energy range 
coincides with the number of electrons Me per lattice site, as for the case of electron-solid phases 
in the conventional two-dimensional electron gas. Analyzing the local density of states at the peak 
energies, we find particular scaling properties of the density patterns if one fixes the ratio vn /Me 
between the filling factor vn of the last partially filled Landau level and the number of electrons 
per bubble. Although the total density profile depends explicitly on Me, the local density of states 
of the lowest peaks turns out to be identical regardless the number of electrons Me. Whereas these 
electron-solid phases are reminiscent of those expected in the conventional two-dimensional electron 
gas in GaAs heterostructures in the quantum Hall regime, the local density of states and the scaling 
relations we highlight in this paper may be, in graphene, directly measured by spectroscopic means, 
such as e.g. scanning tunneling microscopy. 

PACS numbers: 73.43.-f, 73.20. Qt, 73.21.-b, 68.37.-d 

Keywords: Wigner crystal; electron-bubble crystal; high magnetic field; local density of states; scanning 
tunneling spectroscopy 



INTRODUCTION 



As was shown by Wigner in 1934,^ the degenerate 
Fermi gas is unstable towards the formation of a pe- 
riodic triangular lattice of localized electrons (electron 
crystal), once the Coulomb energy prevails over the ki- 
netic one. However, the critical electron density at which 
the transition to the Wigner crystal (WC) occurs is too 
low for usual metals. Nevertheless, the situation is much 
improved if one applies a strong perpendicular magnetic 
field to a two-dimensional (2D) electron gas (2DEG). In 
this case, the single-particle continuous energy spectrum 
is quantized into a sequence of hugely degenerate Lan- 
dau levels (LL's). If one restricts oneself to the electrons 
within the last partially-filled LL, one finds that their ki- 
netic energy is quenched, and the only energy scale is the 
Coulomb energy, which favors the formation of an elec- 
tron crystal at small filling factors ^^i^*^ 

A quantum electron crystal in the presence of a 
disorder potential is expected to become collectively 
pinned and to manifest itself as an insulator.^ While 
at small filling factors the 2D WC with triangular lat- 
tice symmetry^'— is expected to yield the global en- 
ergy minimum, it was predicted that the phase dia- 
gram of the 2DEG includes also electron-bubble crys- 
tals (a periodic lattice with more than one electron per 
site), stripes,— i^i^ and even more exotic quantum Hall 
liquid-crystal phases.^ Unlike electron-crystal phases, the 
prominent quantum liquids, which display the fractional 
quantum Hall effect in the two lowest LLs^^ are transla- 
tionally and rotationally invariant and remain conducting 



even in the presence of disorder. Therefore, these differ- 
ent quantum phases may be distinguished experimentally 
with respect to the behavior in transport measurements. 
For instance, a succession of insulating and conducting 
phases yields a re-entrant integer quantum Hall effect 
(IQHE) in the first^^ (7V=1) and second^^ (TV = 2) ex- 
cited LLs and has been interpreted in terms of a competi- 
tion of such electron-solid and quantum-liquid phases.— 
Further evidence for electron crystals in LLs stems from 
radio-frequency spectroscopy^^ and transport measure- 
ments under microwave irradiation ^^^i^^i^^i^^i^^ which ex- 
cites the collective pinning mode of the electron crystals. 

While all these experimental techniques have been very 
successful in discovering new insulating phases and have 
confirmed a number of theoretical predictions, they are 
indirect evidence for high-field electron crystals based 
on transport measurements - the 2DEG in GaAs het- 
erostructures is buried deep inside the substrate, which 
renders impossible a direct optical observation of a pe- 
riodic electron-crystal lattice, e.g., by means of scanning 
tunneling microscopy (STM).^^ In contrast to the con- 
ventional 2DEG, such optical studies might become pos- 
sible in graphene, a one-atom thick sheet of graphite, 
with unique electronic and mechanical properties In- 
deed, graphene may be viewed as a particular 2DEG, 
where the electrons behave as if they were massless par- 
ticles described by the relativistic 2D Dirac equation. 
In this Dirac equation the Fermi velocity vp plays the 
role of the speed of light c, although it is roughly 300 
times smaller than the latter in vacuum. In addition, the 
Brillouin zone of graphene has two non-equivalent corner 
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points (called Dirac points) which yield a twofold val- 
ley degeneracy and which may formally be described in 
terms of an SU(2) pseudospin degree of freedom. 

In strong magnetic fields, the energy of Dirac fermions 
in graphene is quantized into LL's, the structure of which 
is different from that of non-relativistic electrons in a 
conventional 2DEG. Apart from their unconventional 
(square root) magnetic-field dependence, there exists a 
LL at exactly zero energy, and each LL in the conduction 
band has a counterpart in the valence band. This partic- 
ular LL structure leads to the anomalous integral quan- 
tum Hall effect observed in graphene ^^^^ In addition, 
the energy gap between the subsequent LL's in graphene 
is so large, that it is possible to observe the IQHE even 
at room temperatures.^^ 

As the mobility of graphene samples is further im- 
proved, one may expect to observe the fractional quan- 
tum Hall effect, which has been studied theoretically by 
several authors,^^'^'''^^'^ and also collectively-pinned 
insulating phases such as Winner-crystal and bubble 
phases, as predicted in Refs. ^^'s^illlilll. In con- 
trast to GaAs heterostructures, these electronic phases 
occur at the surface of the graphene sheet and are, thus, 
directly accessible by spectroscopic means. Indeed STM 
has been applied successfully to probe the density distri- 
bution in exfoliated^ and epitaxial^ graphene, as well 
as in graphene on a graphite substrate in a strong mag- 
netic field. This exciting prospect motivated us to cal- 
culate theoretically physical quantities of a 2D electron 
crystal which might be measured in an STM experiment: 
the (integrated) density of states (DOS) and the local 
density of states (LDOS). We should note that the quan- 
tum Hall regime is the only case for which one may ex- 
pect the formation of electron-crystal phases in graphene. 
Indeed, it is predicted that the 2D Wigner crystalliza- 
tion is completely absent in graphene for any electron 
density in the absence of a magnetic field, due to the 
scale invariance of the dimensionless interaction parame- 
ter Ts = /hevp :^ 2/e for a 2D system with a linear dis- 
persion relation. Here, e is the dielectric constant which 
depends on the environment where the graphene sheet is 
embedded. 

In this paper, we discuss the DOS and the LDOS 
for several electron crystals in the N = 2 LL within a 
Hartree-Fock approximation. We have performed similar 
calculations for N = 1, 3 and 4, but we concentrate in 
the present paper on TV = 2 for two reasons. First, the 
DOS and LDOS results for = 2 are representative of 
high-field electron solids - our calculations yield indeed 
similar results for the other LL's. Second, for higher LL's 
there have been no clear indications so far for electron- 
crystal phases in GaAs in the quantum Hall regime. For 
our numerical calculations, we have adopted the itera- 
tive scheme proposed by Cote and MacDonald^*^ which 
has also been applied to calculate the energies and the 
real-space profiles of various electron-crystal phases in 
graphene. As a test of the validity of our code, we have 
corroborated the results obtained in Ref. [3l| and then 



applied it for the calculation of the DOS and the LDOS. 
Notice that, despite the huge amount of Hartree-Fock 
studies of quantum Hall electron-crystal phases, none is 
devoted to study the LDOS of these phases. We have 
calculated the LDOS at energies where the integrated 
DOS has well-pronounced peaks, which fall into two dis- 
tinct classes: bound states at negative energy with re- 
spect to the chemical potential and high-energy peaks 
above. The number of negative-energy peaks is identical 
with the number of electrons Mq per bubble, in agree- 
ment with bubble crystals in the conventional 2DEG;^ 
Furthermore, we find that the sum of the LDOS at these 
Me negative-energy peaks reproduces the real-space den- 
sity profile of the Me-electron bubble crystal. 

This paper is organized as follows. In Section [III we 
outline the basic steps of the Hartree-Fock approximation 
to the 2DEG in graphene. In Section Hill we present nu- 
merical results for the DOS and the LDOS in graphene. 
Finally, we draw our conclusions in Section [iVl 

II. HARTREE-FOCK HAMILTONIAN 

For a partially filled LL N, the low-energy electronic 
properties are captured within a model that takes into 
account states only within this level. In this case, the 
single-particle kinetic energy is the same for all of states, 
and thus only the interaction term is relevant. Further- 
more, we omit the physical spin which we consider to be 
completely polarized, e.g. due to a sufficiently large Zee- 
man effect. The derivation of the Hartree-Fock Hamil- 
tonian for the 2DEG in GaAs has been extensively dis- 
cussed in the literature.— In graphene, the interaction 
Hamiltonian for the 2DEG is similar to that in GaAs, al- 
beit with different form factors due to the spinorial form 
of the wave functions ^^^^ This similarity allows one to 
use the same theoretical methods which were used pre- 
viously to study the 2DEG in GaAs, with the important 
difference that we need to take into account the twofold 
valley degeneracy in the form of an SU(2) pseudospin 
degree of freedom, P = ±1. Provided that inter-LL tran- 
sitions are neglected, we may write the interaction part 
of the full Hamiltonian for the 2DEG of spinless electrons 
in graphene as^^ 

'^int = A[^iv(q)]^P/3,/3(-q)p/3',/3'(q), (1) 

^ /3,/3',q I'll 

where Vc = /Ib^ is the Coulomb energy scale, with 
Ib = \/h/eB the magnetic length, B the magnetic field, 
and e the dielectric susceptibility of the medium, and 
q = {qx^Qy) is a 2D wavevector. The (guiding center) 
density operator in the Landau gauge reads 

X ^ ' 
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Here, A^^ = S/2itI\ measures the LL degeneracy, with 
the square area S of the 2DEG sample, cx,/3 and ^ are 
the electron's destruction and creation operators, respec- 
tively, where X denotes single-particle quantum states 
within the A/'-th LL. Finally, in Eq. ([Tj) the graphene 
form factor J^Ar(q) reads^*^ 



^iv(q) 




^- = 0, 
(3) 



where q= |q|, and Ln{x) is the Laguerre polynomial of 
order n. We note that the 2DEG form factor in GaAs is 
given by^ 



Fn{ci) 



(4) 



It is apparent from Eqs. (|3]) and (|4]) that the graphene 
(relativistic 2DEG) form factor is simply a linear com- 
bination of form factors for adjacent LL's of the non- 
relativistic 2DEG in GaAs. This peculiar fact results 
from mixing the Dirac particle wavefunctions between 
the sites of two sublattices in graphene, and is also a 
consequence of the spinorial nature of these wavefunc- 
tions. Apart from the difference in form factors given by 
Eqs. ^ and (|4]), the 2DEG in GaAs and graphene is de- 
scribed equivalently, as follows from the same analytical 
structure of the Coulomb interaction term given by Eq. 

Finally, we note that the Hamiltonian in Eq. ([T]) is 
SU(2)-invariant with respect to the valley pseudospin. 
In contrast to the physical electron spin, this SU(2) sym- 
metry is approximate. However, SU(2)-symmetry break- 
ing terms are suppressed linearly in q/Ib <C 1 where 
a = 0.14nm is t he ca rbon-carbon distance in graphene 
and Ib = 26/y^[T]nm, i.e. at an energy scale that 
is well below the disorder broadening of the LL's*^*^ 
This physical model is similar to another two-component 
quantum Hall system - if one replaces in Eq. ([1]) ^Ar(q) 
by the non-relativistic form-factor FAr(q), one obtains 
the Hamiltonian for the non-relativistic 2DEG including 
the electrons' spin in the absence of a polarizing Zeeman 
effect. Alternatively, this model may describe a quantum 
Hall bilayer in the theoretical limit of zero layer separa- 
tion, where the two "spin" orientations denote the two 
different layers One may further simplify the model in 
Eq. ([T]) by omitting the valley pseudospin degree of free- 
dom, in which case one presupposes a complete valley 
polarization of the electronic phases, which would maxi- 
mally profit from the exchange interaction. This effective 
U(l) model is described by the interaction term 



(5) 



where the density operator of spinless electrons p(q) is 
obtained from Eq. ([2j) by neglecting the pseudospin in- 
dices. This simplified U(l) model of fully valley-polarized 



graphene which is described by Eq. ([5|) is called U(l)- 
graphene in the remainder of the paper. Now, if one 
substitutes into Eq. ([5]) the non-relativistic form-factor 
Fn{^), one obtains the usual single-layer quantum Hall 
2DEG for spin-polarized electrons in GaAs. 

The Hartree-Fock approximation applied to the 
graphene interaction term in Eq. ([T]) yields^*^ 



(HF) 



I3,Q 

-X^'^{Q)p^^^{Q)], (6) 



where ^ = —(3, and Q's are the reciprocal wavevectors of 
the WC lattice. The Hartree and Fock effective interac- 
tion potentials read, respectively. 



Hid) =- 



Q 

'(Q)=/ ' 
Jo 



l^iv(Q)|V(-Q)(l-^Q,o), (7) 



dxe-^\J^N{Q)\^MxQ)p0,0'{-Q), 



(8) 



where Q = |Q|, Jo is a Bessel function, and the density 
averages are p/3,/3'(Q) = (p/3,/3'(Q)), p(Q) = E/3 P/3,/3(Q). 
We assume a triangular electron lattice for the broken- 
symmetry state, with reciprocal lattice vectors given by 



Q = go 



n n 
2' 2 



m 



(9) 



Here Qo is the length of the basis vector of the reciprocal 
lattice, 



(10) 



UN is the filling factor of the last partially filled LL, 
and Me is the number of electrons per site (Me = 1 
corresponds to the WC, and Me > 2 to an electron- 
bubble crystal with Me electrons per bubble) . The single- 
particle Green's function in the imaginary-time Matsub- 
ara formalisn^i reads 



<^/3i,/32(Q,^^n) 



-N2 



X 



h/kuT 



dr exj[){icjnT) 



lie 



X (^rCx-(|Q,AM4,/3.(0)), (11) 



where T is the temperature, kB is the Boltzmann 
constant, %- denotes imaginary-time ordering, and 
cOn = 7r(2n + l)kBT/h are the Matsubara frequencies. 
^/3i,/32 (Q7 ^^n) may be determined self-consistently from 
the quadratic Hamiltonian (|6]) by using the Heisenberg 
equations of motion within the iterative-solution method 
proposed in Ref. UBI which we adopt in the present work. 
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After analytic continuation to real frequencies icOn 
cj + iO^, G/3^,/32(Q, icc;^) yields the retarded Green's func- 
tion which may be used to calculate the DOS g{uj), 

g{ij) = -N-^iT Im Gf3^f3{Q = 0, ic^n ^ ^ + ^0+), 

/3 

(12) 

and the LDOS A(r,cj), 
A(r,cj) = -N-^^lm G/3,/3(r, icj^ ^ cc; + zO+), (13) 

where the Green's function in real space reads 

G'/3i,/32(r,icc;n) = (27r/|)"^ ^exp(-iQ • r)J^Ar(-Q) 

Q,/3 

X G/3i,/32(Q,^^n). (14) 

III. RESULTS AND DISCUSSIONS 

In this section, we discuss the spectroscopic properties 
of the electron-solid phases for the LL = 2. As al- 
ready mentioned in the introduction, we concentrate on 
this LL for illustration purposes and because they are 
mostly significant from the physical point of view. We 
have obtained similar results for A" = 1, 3 and 4 (not 
discussed here). 

We have chosen three different electron-solid lattices 
which have the same ratio un/Mq = 0.14 ^ 1/7, and 
hence the same lattice period given by Eq. (p!Q|) . Our 
choice for the vn /Mq ratio is rather arbitrary. We note 
that the Me = 1 and Me = 2 states yield in graphene the 
global energy minima, while the Me = 3 does not. It has 
been shown in Ref. [3l| that the ground state of graphene 
at vn < 0.43 is an anisotropic Wigner crystal whereas at 
0.28 <vn < 0.43 the ground state is the Me = 2 bubble 
crystal, and at vn < 0.28 the Wigner crystal yields the 
lowest energy (Me = 1). Therefore, the Me = 3 phase 
is not the lowest-energy state in graphene; nevertheless, 
it is useful to analyze on the same footing all three cases 
Me = 1,2,3. 

For completeness, we mention that we have also cal- 
culated the cohesive energies of other types of electron- 
crystal phases (not only triangular bubble phases, but 
also anisotropic Wigner crystals). Our results for the 
energies coincide with those of Ref. [3l| with excellent ac- 
curacy and, therefore, corroborate the DOS and LDOS 
results discussed below. 



A. (Integrated) Density of states 

Our results for the DOS in graphene at A" = 2 are 
presented in Fig. [TJ 

We find that the DOS consists of two well-separated 
classes of peaks: well-defined low-energy peaks are found 
below the chemical potential /i, which is shifted to zero 



energy, whereas the large number of peaks above /i are 
not that easily distinguished. We note here that in 
graphene the number of low-energy peaks in all cases 
is equal to Me, the number of electrons in a bubble. 
The same result has been obtained before in Hartree- 
Fock studies of the simpler single-layer 2D quantum 
Hall system.^ We checked that the same property holds 
true also for U(l)-graphene, and in non-relativistic two- 
component quantum Hall systems, such as a bilayer with 
zero layer separation. 

In the simpler single-layer 2D quantum Hall system 
in GaAs, the DOS at Me = 1 exhibits the features of 
the Hofstadter butterfiy structure.^ It means the follow- 
ing: given that the filling factor may be represented by 
a ratio of two integers p, q without a common divisor, 
= these integers p, q determine then the struc- 
ture of the single-particle energy spectrum of the system; 
namely, there should exist p low-energy levels and q — p 
high-energy levels (Hofstadter butterfiy counting rule). 
In the DOS, which is a function of frequency, these en- 
ergy levels are recognized as smoothed peaks. The Hofs- 
tadter butterfiy counting rule was confirmed for Me = 1 
in the single-layer 2DEG in GaAs, while for Me > 2 it is 
claimed that counting of the single-particle levels is dif- 
ferent: the number of low-energy peaks is equal to Me, 
whereas nothing is known about what is the precise rule 
for counting the number of high-energy peaks. ^ 

(A) (B) (C) 

Bound states 
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Single-particle excitations 
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FIG. 1: Logarithmic plots for the density of states g{ijj) of 
graphene in the A = 2 Landau level at (A) Me = 1, z^at = 
0.14; (B) Me = 2, i/iv = 0.28; (C) Me = 3, z/jv = 0.42. 
The frequency uo is given in units of the Coulomb scale Vc- 
The a; = frequency is the position of the chemical potential 
(Fermi energy) /i. The infinitesimal imaginary frequency shift 
uj cj + iO^ is approximated by a; ^ cj + i^o;, with 5uj — 10~^. 
Figures (A)-(C) in the first row yield the DOS in the low- 
energy frequency range (bound states of electrons), while in 
the second row the high-energy peaks of the DOS correspond 
to single-electron excitations above the ground state of the 
lattice. 

What will be important in the following discussion is 



the order of indexing of the DOS peaks. We will count 
the peaks in the DOS with respect to increasing the fre- 
quency u. In Fig. mA), the first DOS peak is obvi- 
ously the lowest-energy one with energy ^ — 0.24Fc- The 
second peak in (A) is a higher-energy one with energy 
^ O.llVc' The numbering of peaks continues until we 
reach the utmost-right peak with energy ^ 0.3Vc. The 
same indexing rule is applied to the cases (B) and (C). 
One should note that in all three cases the DOS peaks 
with the same index may have rather different energies: 
while in (A) the second DOS peak belongs already to the 
high-energy region, in (B) it still lies below the Fermi 
level. 

In addition, one should have a procedure of extract- 
ing the energies of the DOS peaks from the smoothed 
DOS vs frequency dependence shown in Fig. [H In the 
low-energy regime it may be always done reliably. In the 
high-energy regime, however, there is a larger number 
of closely-located DOS peaks, the shapes, widths, and 
amplitudes of which depend sensitively on the imaginary 
frequency shift 5^. The latter is used for the analyti- 
cal continuation into the upper complex half-plane of the 
Green's function, iuon ^ a; + i5^^. Physically, this imag- 
inary frequency shift represents a level broadening due, 
e.g., to disorder. We have found that the best way to ex- 
tract only those peaks which are physical is to place a cut- 
off A on the DOS peak amplitude, so that peaks with am- 
plitude less than ^maxA are neglected, with ^max ^ 
the maximum peak amplitude. In our study, we have 
chosen 5^ = 10~^, and A = 0.5. The number of shells of 
reciprocal lattice vectors Q's is A/'sh = 8, so that the ac- 
tual number of vectors is Nq = 241. Single-particle ener- 
gies which are extracted from the smoothed DOS, will be 
used below in the calculation of the LDOS. In the U(l)- 
graphene and the single-layer cases at the same densities 
considered here, we are able to extract 7 DOS peaks; at 
UN = 0.14, Me = 1, there is one lowest-energy peak, and 
6 high-energy ones. These are exactly the numbers of 
single-particle levels dictated by the Hofstadter butterfly 
counting rule.^ 

In graphene and the quantum Hall bilayer, studied 
within the two-component model, we obtain the num- 
ber of identified DOS peaks around 14 [with deviation 
of not more than one wrongly identified peak]. Due to 
the additional SU(2)-symmetry in the latter two cases, 
it is natural to expect that the number of single-particle 
levels is thus doubled. 



B. Real-space density profile 

For later comparison with our results for the LDOS, 
we calculated the real-space electron density profile 

n{r) = ^^exp(-zQ.r)^^(-Q)/)(Q) (15) 
^ ^ Q 

for the same choices of {vn^ Mq) as in Fig.[Tl The results, 
which are shown in Fig. [21 agree with previous calcula- 
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tions for graphene performed by Zhang and Joglekar*^ 



(A) (B) (C) 

[0.0,2.1] [0.3,1.7] [0.6,2.4] 




FIG. 2: (Color online) Real-space density profile n(r) in 
graphene in the N — 2 LL. Choices (A)-(C) are the same 
as in Fig. [T] Minima and maxima of n(r) written inside the 
square brackets as [min, max] correspond to the values [0.0, 
1.0] in the colour plots (blue and white colours, correspond- 
ingly). 



C. Local density of states 

Our results for the LDOS in graphene are presented 
in Figs. [31 H and [5l The LDOS patterns is plotted for 
all three cases (A)-(C), as in Fig. [U and at the ener- 
gies of all extracted single-particle DOS peaks situated 
in increasing order. We obtain that the rescaled (to the 
range of [0.0, 1.0]) real-space patterns of A(r,cc;), calcu- 
lated at the first four DOS peaks for all three choices 
of {un-, Me), coincide among themselves. There is also 
an approximate mapping between the LDOS patterns at 
the fifth DOS peak, although less pronounced than for 
the first four ones (one sees correspondence between the 
positions of maxima and minima, but the colors deviate 
slightly in each case). For larger values of the peak in- 
dex, we start to see considerable discrepancies between 
the LDOS patterns. Also the number of extracted peaks 
A^p is different for each case. The latter property is due 
to the very approximate nature of our extraction proce- 
dure: while low-energy peaks are always identified reli- 
ably, the high-energy peaks are determined only approx- 
imately. However, the accuracy is quite good. We also 
note a very interesting property of the LDOS at the last 
two peaks for (B)-(C): the LDOS patterns are identical, 
but their positions are swapped. We do not have any 
physical argument why this should be the case, but it 
could be a feature that appears when M > N. However, 
this statement is a mere speculation and a more detailed 
investigation is required to clarify this aspect. 



D. Comparison with the real-space density 

Now, to compare the LDOS patterns shown in Figs. 
[31 m [5] with the real-space density profile n(r) defined in 
Eq. ([15]) and plotted in Fig. [21 we introduce the resummed 
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FIG. 3: (Color online) LDOS A{r,Lj) for graphene at un = 
0.14, Me = 1 [case (A)]. Contour colours are graded in the 
same way as defined in Fig. O The contour plots are ordered 
with respect to the index of DOS peaks [indicated above the 
plots]. The number of extracted DOS peaks is A^p=13. 



FIG. 5: (Color online) LDOS A(r,cj) for graphene at un 
0.42, Me = 3 [case (C)]. A^p=14. 




FIG. 4: (Color online) LDOS A{r,oo) for graphene at un 
0.28, Me = 2 [case (B)]. Np = 15. 



LDOS A(r,cj), defined for a fixed single-particle energy 
uji as a sum of all LDOS patterns at smaller peak energies, 




FIG. 6: (Color online) Resummed LDOS A{r, cu) for graphene 
at the three first DOS peaks for (A) [un = 0.14, Me = 1]. 
Contour colours are graded in the same way as defined in Fig. 

m 



where the sign ^ means mapping between the rescaled 
to the [0.0,1.0] interval quantities. This means that the 
real-space density of the Me-electron bubble crystal is 
determined by the sum of the LDOS at the Mq negative- 
energy peaks. More surprisingly, because of the corre- 
spondence between the LDOS patterns of the low-energy 
peaks for all different Me bubble crystals, one may de- 
termine the real-space density pattern of the Me bubble 
crystal by summing the LDOS of the Me peaks of lowest 
energy for amy of the electron-solid phases - the LDOS 
patterns of the Me = 1 Wigner crystal, e.g., contains 
thus the information of the density of all other Me bub- 
ble crystals. 



(16) 



Given the excellent coincidence of the LDOS patterns 
shown in Fig. [6] with the real-space densities in Fig. O 
one may empirically write 

n(r. Me = ^ i(r, cj^), i = 1, 2, 3, (17) 



E. U(l)-graphene: local density of states 

In Fig. [71 we present the LDOS for U(l)-graphene. The 
number of extracted DOS peaks for all three density 
choices (A)-(C) is Np = 7. This is in accordance with 
the Hofstadter butterfiy counting rule. We also see an 
excellent correspondence of the LDOS for the first four 
DOS peaks, then, also a good coincidence of the LDOS at 
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the peaks 6-7 for (A)-(B), whereas these two LDOS pat- 
terns for (C) interchange places, as compared with the 
patterns 6-7 for (A)-(B). This interchange phenomenon 
is the same as observed for graphene and is not yet un- 
derstood. 

In general, the U(l)-graphene results coincide numeri- 
cally with those for graphene when one takes into account 
the SU(2) symmetry for the valley pseudospin. This in- 
dicates that one may use the U(l)-model instead of the 
more complex SU(2)-symmetric one for the discussion of 
the density patterns, hence simplifying further calcula- 
tions on graphene. Moreover, it indicates that in the 
electron-crystal phases considered above the valley de- 
gree of freedom is fully polarized. 




FIG. 7: (Color online) LDOS A{r,oo) for U(l)-graphene. The 
number of extracted DOS peaks A^p = 7. 



graphene is given by the number of electrons per site Mq. 
This result is similar to the previous DOS calculation in 
the Hartree-Fock approximation for GaAs.^ 

We found that the rescaled LDOS is identical for differ- 
ent filling factors un, as long as the ratio un/Mq, which 
determines the lattice spacing of the Me-electron bubble 
crystal, is kept fixed, and the LDOS frequency is taken 
at the DOS peak with the same index (for the first four 
indices). In particular, this result yields an unexpected 
conclusion that, e.g. by fixing the filling factor = 0.14 
in the TV = 2 LL, and using STM, one could observe in 
the LDOS the whole succession of electron-crystal den- 
sity patterns with Me = 1, 2, 3 by fixing the applied 
STM voltage at the consecutive first three single-particle 
excitation energies, and summing up the LDOS to obtain 
the resummed LDOS A{r^u;). 

We believe that this LDOS correspondence holds true 
for all single-particle excitations resolved as individual 
DOS peaks so far (accounting for interchanging of the 
last two peaks in the Mq = 2,3 cases) and for all LL's 
(similar conclusions follow from our calculations in the 
LLs N = 1, 3, 4). 

We also obtained the same LDOS correspondence for 
other models of the 2DEG: (i) in a single-layer GaAs 
heterostructure; (ii) U(l)-graphene; (iii) bilayer. This 
implies that the observed LDOS un/Mq scaling is inde- 
pendent of the underlying interaction potential and the 
number of inner discrete degrees of freedom. The fact 
that the U(l)-graphene results coincide numerically with 
those for graphene indicates that the electron crystals 
considered here are completely valley-pseudospin polar- 
ized. 



IV. CONCLUSIONS 

The aim of this work is to show how a high-field 
electron-solid phase in the 2DEG may be detected by 
optical means in graphene. We have calculated the DOS 
and the LDOS of electron-solid phases in the Hartree- 
Fock approximation in the N = 2 LL. 

We show that the number of low-energy DOS peaks in 
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